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Abstract. Coupling between heat and electrical currents is at the heart of thermoelectric 
' processes. In a thermoelectric system this may be seen, from a thermal viewpoint, as an 

S additional thermal flux linked to the appearance of an electrical current. Since this additional 
I ■ flux is associated with the global displacement of charge carriers in the system, it can be qualified 

' as convective in opposition to the conductive part related to both phonon transport and heat 

, transport by electrons under open circuit condition as, e.g., in the Wiedemann- Franz relation. 

O ' In this article we demonstrate that considering the convective part of the thermal flux allows 

O \ both new insight into the thermoelectric energy conversion and the derivation of the maximum 

power condition for generators with realistic thermal coupling. 

> 

. 1. Introduction 

The field of thermoelectricity was born in the 1820s after Seebeck reported his observation 
^ ■ of a phenomenon he interpreted as thermomagnetism. Orsted, in the same years, gave a 

^ ! better interpretation of the phenomenon and termed it thermoelectricity: the appearance of 

an electrical current in a circuit subjected to a temperature difference. The name of Seebeck 
was later given to the coefficient linking the temperature difference AT and the resulting 
electromotive force AV^: a = —AV/ AT. A few years after Seebeck's works were published [T], 
Peltier noticed that, conversely, an electrical current / may absorb or release heat at a junction 
. between two different materials [2j. The flux of absorbed heat Iq is related to the electrical 

current / through the difference of Peltier coefficients U between the materials of the junction: 
Iq = (n^ — Ila)I (with / positive when flowing from material a to material b). The relation 
between these two effects was later clearly expressed by Thomson who established the connection 
between the Peltier and Seebeck coefficients: 11 = q:T, with T the temperature of the considered 
material [3J. This relation demonstrates that both effects have the same physical origin: the 
ability of charge carriers to carry heat along with electrical charge. As thermoelectric (TE) 
devices possess the particularity to strongly couple heat and electrical transports, they can be 
used in refrigeration or energy conversion applications [4J. However to achieve good performances 
both TE materials and device design have to be optimized. For materials, optimization amounts 
to increasing the so-called figure of merit ZT (see for example the recent review by Shakouri 



[5]). Yet, to make the best use of the materials it is necessary also to lead a reflection on the 
integration of the devices into an external environment. One thus needs a complete model of 
thermoelectric device to fully comprehend the interaction between the constitutive laws of the 
device and the laws related to the external exchanges. 

In this work we propose a simple model of thermoelectric generator (TEG) that accounts 
for the heat transported by electrical carriers; this results in the introduction of an effective 
thermal conductance associated with this phenomenon. As this heat propagation is associated 
with global movement of particles, by analogy with fluid mechanics, we call it convection. We 
then study and illustrate the consequences that this additional term may bear on both the 
physics of TE energy conversion and on the practical design and optimization of thermoelectric 
generators (TEGs). 



2. Thermoelectric modeling: Accounting for a convective term 

Let us consider a thermoelectric generator connected to two heat baths at temperature T^ot 
and Tcoid (^hot > T'coid) respectively with an average temperature T = (Thot + 2coid)/2. 
Though a thermoelectric device may operate indifferently either as a cooler or as a generator, 
in this work, we restrict our analysis to the case of a generator to which an electrical load 
resistance Rioad is connected in order to extract power. Now that we have set the thermal and 
electrical environment, we state the phenomenological laws that govern the device by relating the 
electrical current / and the average thermal current Iq to the voltage AV and the temperature 
difference AT across the generator. The equations follow directly from the force- flux formalism 
introduced by Onsager and then extended to thermoelectricity by Callen 16]. For our purpose, 
we find convenient to use the macroscopic coefficients and potential differences, rather than local 
coefficients and gradients, and we get: 

I=^AV + ^AT (1) 
aT ( c?T \ 

Iq = —av^\—-^kAat (2) 

where R is the electrical resistance, Kg the thermal conductance under open-circuit condition 
and a the Seebeck coefficient. The TEG is also characterized by its figure of merit ZT obtained 
from these three parameters: ZT — a^T/RK^). This dimensionless number gives a direct 
quantitative estimation of the generator's quality p^J. The convention for the currents and 
potential differences, AT^ — Vi — V^ and AT = T^ot — ^coid? ^ire displayed in Fig. ([1]). 

Equation ([1]) gives a definition of the Seebeck coefficient under open-circuit condition (/ = 0), 
namely a — — Ay/AT. If the Seebeck coefficient a vanishes, thermoelectric effects become 
negligible, and one recovers the traditional macroscopic expressions of Ohm's law and Fourier's 
law, from Eqs. ([T]) and (E]) respectively. 

The two previous equations may be combined to express the thermal flux as the sum of two 
terms: 

Iq = aTI + AT (3) 

The first term is the contribution due to the heat transported within the electrical current, 
the so-called convective thermal flux [HI |9] , and the second is the conductive part traditionally 
associated with Fourier's law. The quantity of energy transported by each carrier is |a|Te, e 
being the elementary electric charge. 

To compare these two terms we now rewrite Eq. ([3]) to obtain a form similar to Fourier's law: 
A relation of proportionality between the heat current and the temperature difference. For this 
purpose we express the electrical current through the resistive load as a function of AT. From 
Ohm's law: / — AlZ/i^ioad? we get: 




Figure 1. Thermoelectric generator: Decomposition of the TEG thermal conductance into a 
conductive part Kq and a convective part i^conv 
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After substitution of the above expression for / in Eq. ([3]), the average heat flux may be expressed 



as: 



\ 



-Rload + R 



AT = KtegAT, 



(5) 



where an additional term, A^conv, foi" the total thermal conductance of the generator now appears 
clearly. This term is directly related to the electrical condition on the circuit since the load 
resistance i?ioad is explicitly involved in its definition. The effective thermal conductance of 
the generator K^eg is thereby equivalent to the two thermal conductances A^conv and Kq 
in parallel as illustrated in Fig. ([1]). When i?ioad ~^ co, i.e., in open-circuit condition, we 
recover, as expected, the fact that K^eg is equivalent to Kq only, whereas when i?ioad ~^ 0, 
KTEG-a^T/R + Ko. 

To derive Eq. ([5j) we made a few assumptions: We considered only an average value for the 
heat flux; we neglected the Joule effect and the thermal energy converted into electric energy 
along the device. These approximations are reasonable if the temperature difference AT is kept 
small; a condition that must be satisfied to make use of constant coefficients to describe the TE 
generator. All these assumptions amount to stating that we only consider a linear framework. 
Second, we neglected the heat flowing through the electrical load. So, dealing with a true unileg 
generator as in Ref. [TU], implies A^teg ^ ^load- This configuration may also be related to a 
generator with two legs: The parameters then used are only equivalent ones [7J, in which case 
the load sustains no temperature difference, so that no heat flux goes through it. 



3. Illustrative examples of the effects of thermoelectric convection 

In this section we illustrate the impact that the convective part of the generator thermal 
conductance may have on the properties and operation of a TEG. We first look at the 
construction and meaning of the figure of merit ZT; next we discuss the question of optimization 
of thermoelectric generators when the connections to the heat reservoirs are not perfect; finally 
we examine the effective thermal conductance of two thermoelectric modules in parallel. 



3.1. On the meaning of the figure of merit ZT 

To discuss the properties of thermoelectric systems, we find convenient to proceed by analogy 
using the classical heat transfer description in terms of convection and conduction. To determine 
whether heat transfer in a fluid is dominated either by conduction or by convection, one may 
introduce the Prandtl number cTp defined as 

- ^ (6) 

where v is the kinematic viscosity of this fluid, characterizing momentum diffusion, and D 
its thermal diffusivity, characterizing heat diffusion [llj. We retained this idea to define a 
thermoelectric Prandtl number (crp)^^ given by the ratio of the convective thermal flux and of 
the conductive thermal flux: 

^ KconvAT ^ a'T ^ ZT 

^^P^™ KoAT Ko(i? + i?load) l + i?load/i? 

It appears that (crp)^^ is proportional ZT. The proportionality coefficient depends on the ratio 
between the electrical resistances of the generator and the load. If the electrical circuit is shorted, 
i.e., if i?ioad <^ Ri the thermoelectric Prandtl number (crp)^^ is then equal to ZT, which is the 
maximum value of the ratio between the convective part and the conductive part of the thermal 
flux. It is consistent with the fact that ZT must be increased to improve the TEG performances: 
The convective thermal flux is a direct consequence of the coupling between heat and electrical 
charge transport, and, as such, it is a necessary process to obtain energy conversion. Conversely, 
the conductive thermal flux amounts to pure loss effects only: So maximizing ZT is equivalent 
to maximizing the useful part of the thermal flux with respect to the losses. 

The ratio i?ioad/^ controls the convective component of the heat flux. However, maximization 
of (crp)^g does not lead to the best performances: To extract power from the heat that flows 
in the generator, one has to satisfy impedance matching too. If i?ioad becomes too small, 
the electrical current, and hence the convective thermal current, will increase but the power 
will decrease. It appears that the optimal value of (crp)^^ associated to the maximum power 
output is not ZT but ZT/2. We can even extend this analysis to the maximization of the 
TEG efficiency: As shown by loffe [7J, for a TEG with constant parameters the efficiency is 
maximized when i?ioad/^ — \/l + ZT. Thereby, for the efficiency the optimal value of (crp)^^ is 
ZT/(1 + Vl + ^T). 

This last result should be compared with the compatibility approach of thermoelectric 
optimization: Snyder and Ursell [12j demonstrated that in order to obtain the best efficiency, 
the relative current u — J/ (avVT) (with J the electrical current density and n the local thermal 
conductivity under open circuit-condition) has to be kept constant along the device and should 
be equal to an optimal value s — + ZT — l)/aT. Expressed differently the condition for 
efficiency maximization reads: aTJ/{i<iVT) — + ZT — 1). As the term on the left hand 
side is a local definition of (ap)^^, we recover exactly the same result as above. It is interesting 
to note that similar conditions were derived using the ratio of the total heat current to the 
convective thermal flux in Refs. [131 El- 

Finally the notion of convective thermal conductance can also be useful to practically 
determine the figure of merit of a thermoelectric generator: In 1966, Lisker proposed to 
obtain ZT by measuring the effective thermal conductance of the TE module for two different 
electrical load conditions instead of extracting from three different measurements the value of 
each parameter i?, a and [15j. The ratio K/^y=()/ is indeed equal to + 1. 



3.2. Thermal matching of thermoelectric generators 

Assume that the generator is not directly coupled to the heat reservoirs; the finite thermal 
conductance of the contacts must be taken into account. These are denoted A^hot on the hot 
side and A^coid on the cold side of the generator. The total equivalent conductance for both 
contacts is given by i^(^ntact ^ ^hot + ^coid- "^^^ main effect of these imperfect contacts is that 
the actual temperature difference across the TEG is not AT but AT' = T^m — 2cM as illustrated 
on the left side of Fig. 1^. 

We first need to determine the relation between these two temperature differences. To do so 
we use an analogue to the potential divider formula for the temperature, assuming that the heat 
flux remains constant along the structure (this hypothesis is the same as that used to derive the 
expression of A^conv)- We obtain the following relation: 

AT' = ThM - TeM ^ ^ ^7'^^' AT (8) 

^TEG + ^contact 

Since the thermal conductance A^teg depends on the electrical load resistance through A^conv? the 
temperature difference across the generator, AT', is modulated by the electrical load condition. 
The consequence of this modulation is that the electrical circuit modeling of the TEG as diplayed 
in Fig. ([2|) is no longer a Thevenin generator: In such model the voltage source must be 
independent of the load which is not the case when finite thermal contacts are introduced. 
To keep using a proper Thevenin generator we split the thermoelectric voltage Vqc = alST' 
into two terms: One is independent of the electrical load, and the other is linked to an internal 
resistance [16j: 

TT Arri ^contact J- 

Voc = a AT I-^ —j^ (9) 

^0 + ^contact ^7=0 ^contact 

This equation is of the form: Vqc = ~ ^^th = a^T/{K^^Q + ^contact)- We thus 
obtain a rigorous Thevenin modeling of the electrical part of the TEG with the definitions of 
the open circuit voltage given by V^^ and the internal resistance is i?TEG — R + R''. As for 
the effective thermal conductance, the effective electrical resistance of the generator shows an 
additional contribution directly related to the thermal conditions imposed on the generator 
through the thermoelectric coupling coefficient a. Further, it is remarkable that the expressions 
of R' and A^conv present the same form. One can switch from one to the other by changing 
thermal conductances to electrical resistances and vice versa. Recently, the introduction of an 
extra electrical resistance depending on thermal conditions was also suggested by Spry [17j . 
With our model, the power produced by the TEG may thus be simply expressed as: 

P^AVI^ ^oc'-Rload ^^Q^ 
(i?TEG + i?load)2 ^ 

from which we deduce that power maximization is obtained for i?ioad — Rteg- Such a condition 
was already derived mathematically by Freunek and co-workers [18] . 

As regards thermal matching, Stevens proposed that to maximize power production, the 
TEG should be thermally matched with Kq = ^contact- Let us examine the situation using 
the equations we just derived. An expression for the maximum power may be obtained using 
Eq. (fTOl) in the case of electrical matching and the definition of i?TEG- 

^ _ (KcontactAr)2 ZT ^^^^ 



4{K,^, + i^contact)T 1 + + Kcontact/i^z=o 

For a given value of ZT, i^max shows a clear dependence on the ratio i^j^o /^contact • If we 
consider a fixed value of A^contact? due for example to technological constraints, the choice of the 



Figure 2. Thermoelectric generator with dissipative thermal coupling (left) and maximum 
power output as a function of the ratio i^o/^contact foi" different values of ZT (right). 

value for K^^^ will strongly impact the ability of the generator to produce power. On the right 
side of Fig. ([2]) we show the relation between Pmax and A^j^q /^contact foi" three different values 
of ZT (1, 3 and 10). Note that for each point electrical impedance matching is achieved. We 
see that this condition for power maximization through thermal impedance matching slightly 
differs from the one obtained by Stevens, and that it depends on ZT. Actually the optimum 
value of Kq corresponds to the thermal matching between A^contact ^ind A^teg- ^teg — ^contact 
[16j. The dependence on ZT is explained by the fact that the greater ZT is, the greater the 
contribution of A^conv to A^teg is. Thus for sufficiently high TEG efficiency, i.e., high ZT, the 
discrepancy between these two conditions of thermal impedance matching is not negligible: See 
the difference in Fig. ([2|) between the maximum power and the value of power for the condition 
Kj^Q = ^contact displayed by the vertical dot-dashed line. 

Finally we express the simultaneous conditions of thermal and electrical matching in terms 
of ratios i^o /^contact and R\oad/R (as R and Kq are more suitable quantities for device design 
than Kteg and i?TEG)- 



Once again, we emphasize the symmetry that arises between thermal and electrical quantities. 
In a recent publication, Yazawa and Shakouri obtained the same relations from a more 
mathematical viewpoint [20j . 

3.3. Effective thermal conduction increase in inhomogeneous material /^7]/ 

Our last example concerns the determination of the effective thermal conductance of two 
thermoelectric modules in parallel as represented in Fig. We keep the same notations 

for modules' properties except that we add a subscript (1 or 2) to distinguish the two modules. 
First, we derive the total average thermal flux Iq^ inside the system: It is given by the sum of 
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Figure 3. Association of 2 thermoelectric generators in parallel. 



the average thermal flux for each module, with a conductive and a convective part for each. So 
we have: 

^e, = ^1 + ^2 = ^i^^i + i^o,i AT + a2Tl2 + Ko,2AT (14) 

As we focus on the open electrical circuit condition to determine the equivalent thermal 
conductance i^o,eq5 the total electrical current is set to / = 0. However this does not mean that 
the currents Ii and I2 vanish as well: The only consequence is that Ii = —l2- A simple analysis 
of the electrical circuit in Fig. ([3|) leads to the following expression for the internal current: 

This internal current only exists if the electromotive forces, given by a AT, do not compensate 
each other, i.e., if ai 7^ 0^2- In that case we have to consider the convective contribution to the 
thermal flux. Substituting Eq. (|T5]) in Eq. (|T^ and simplifying by AT we get the equivalent 
open circuit thermal conductance: 

i^o,eq = Ko^i + Ko,2 + ^"'""'/^ (16) 

ill -h 1X2 

This last equation demonstrates that the convective component of the thermal flux has a 
significant impact on the effective thermal conductance of the two thermoelectric modules in 
parallel if they possess very dissimilar Seebeck coefficients. However this term is of limited 
interest if the electrical resistances of the modules are too high: the magnitude of electrical 
internal current is then small and, consequently, so is the convective thermal flux. Knowledge 
of the effective thermal conductance may be useful when dealing with superlattices for which 
the layers can be viewed as parallel modules [9J. More surprisingly. Price obtained a similar 
expression for the effective thermal conductance of an ambipolar conductor where the subscripts 
1 and 2 are associated with electrons and holes repectively: The two types of carriers are indeed 
considered as evolving in parallel inside the system [22]. 



4. Conclusion 

Convective thermal flux will have more and more impact on the behavior of thermoelectric 
generator as the performance of materials will be increased. Indeed we have shown that the 
proportion of convective to conductive thermal current is directly related to the figure of merit 
ZT which is expected to increase significantly in the future to allow a wider use of thermoelectric 
power production. In this paper we have demonstrated that this additional contribution to heat 
flow should not be neglected for device design as well as for materials research, especially those 
concerning composites where internal currents can develop and lead to a performance decrease. 
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